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Abstract 

The decay of the Higgs boson into two photons is discussed. Using Sudakov tech- 
nique, we study a new type of double logarithms in the case of the light quark loop. 
First, we examine the origin of such logarithms at the one loop level. Then we 
perform a two-loop analysis, and finally we present an explicit result of the resum- 
mation of double logarithms in all order in QCD. The phenomenological applications 
are discussed. 
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1 Introduction 



The neutral scalar Higgs boson H is an important ingredient of the Standard Model (SM) 
and is the only SM elementary particle which has not been detected so far. The lower 
limit on tuh of approximately 70 GeV was obtained from direct searches at LEP Q. 
The one sided 95% confidence level upper limit on itlh is 420 GeV [I]. The perturbative 
unitarity suggests an upper limit on the Higgs boson mass to be about 1 TeV @. 

The properties of the Higgs particle were studied in details by many authors ( see, for 
example reviews || [|, |, [7|, §, || and references therein) . 

The decay of the Higgs boson into two photons will play an important and unique 
role. It is of interest for the following reasons: 

• for the light Higgs masses this decay mode provides a clear signature for the search 
at hadron colliders 

• the 77 width determines the cross section for the Higgs production in 77 collisions 
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• it is well known that the H — > 77 vertex can serve as a counter of the particles with 
masses larger than the Higgs boson mass. 

It is widely accepted that a detailed experimental investigation will be performed at the 
Next Linear e + e Collider which will also operate in the photon-photon mode |H| , as 
a Photon Linear Collider ( see the review of physical applications in [14j]). The possibility 



of studying scalar Higgs boson in the s-channel production, for example in a reaction 

1 + (1) 

is not so obvious, because of large background from direct quark-antiquark production. 
Fortunately, it was demonstrated that using the polarizations of initial photons one can 
suppress background process at high energy despite the QCD corrections to the back- 



ground process is large |TT], [12], y|. If the task of efficient suppression will be achieved, 
one can hope that the Higgs will be investigated in detail in 77 collisions and the idea of 
precise determination of the #77 vertex will then become an important phenomenological 
issue. 

Therefore, the #77 vertex should be investigated as thoroughly as possible in the SM 
as well in the Supersymmetric extension of the SM (SSM). 

In the present paper we study QCD corrections to the decay rate in the limit of the 
small ratio of the quark and Higgs masses, m q and In this limit, an interesting type 
of QCD double logarithms (DL) appears and we make an attempt of their resummation. 
The fact of the emergency of DL in the process H — > 77 was pointed in [ 16| . The similar 



phenomena was found and studied in |T7|, [15| in the another physical process, 77 — > qq in 



the Jz = channel. The authors of |L5] have performed a careful analysis of the origin 
of the DL appearance and studied DL terms up to two-loop accuracy. 

There are two reasons why the problem of DL in H — > 77 seems to be worth of 
being investigated. The first reason is more or less academical. Once we have found a 
process where a new type of DL (in a sense that it is not a Sudakov-like logarithms) 
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appear it is interesting to understand their origin and even resum them. From QED and 
QCD we know a few examples where DL appear and can be resummed in observable 



quantities (see, for example, reviews [JTSJ, [HJ, and also recent paper [TJj and references 
therein). We claim now that the H — > 77 decay considered here is the next new example. 
The second reason of interest is phenomenological. As we mentioned above the precise 
measurement of the H — > 77 decay width will be performed at Next Linear Collider. 
We note at the very beginning that in the SM the contribution of light quark to the 
amplitude is proportional to the mass of this quark and therefore is small. However 
with an improvement of precision, the contribution of heaviest from 'light quarks' , the 
b-quark , can become visible in future precise measurements at NLC. In addition we want 
to stress that in the SSM the coupling of the b quark with the Higgs may be large due to 
an additional large factor tg{(3) and that the contribution of the b quark at large tg(j3) 
may be as important as the contribution of top quark or W boson. 

In present paper we start with the one-loop analysis and demonstrate that one-loop 
and two-loop contributions contain a new sort of double logarithms 

£ = a s log 2 ' ] 



\™>hJ 

If the ratio m q /mH is quite small the parameter £ can be large enough and the QCD 
expansion in £ is not valid anymore. For example, considering a situation when the Higgs 
has a mass mj? = 500 GeV and decays into two photons through 6-quark loop and b 
quark has a mass m& = 4.5 GeV, we see that indeed £ ~ 1 and resummation of double 
logarithms is necessary. In the present paper we show that this goal becomes feasible 
when the physical origin of DL is understood and QCD series are properly organized. 
First, we examine an origin of such log's by using Sudakov technique and then perform 
a two-loop calculation in DL approximation. The resummation of DL terms can be done 
and we present an exact answer in DLA. 

The paper is organized as follows. In the section 2 we discuss the one-loop ampli- 
tude and the two-loop QCD correction. Section 3 contains the resummation of the DL 
contributions in all orders the perturbative QCD and a brief discussion. 



2 One and two loop results 

We start with considering the amplitude of the process H — > 77 at leading one-loop level. 
The amplitude of the Higgs decay into two photons through quark loop can be presented 
in the following form 

M ( q uarks ){H ^ = (G F V2) 1/2 ^N c ^2Q 2 q F( Xq ), (2) 

where the tensor is defined as 

d>» = (hk 2 )g^ - k$k v u (3) 

ei,e2 and ki,k 2 are polarization vectors and momenta of the photons, N c - number of 
colors, Q q - the electric charge of the quark q and the variable x q is defined by the 
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following relation: 



x n 



m 



(4) 



H 



The corresponding decay width was calculated in |20], ^T[] at the one-loop level. In our 
notations, the quark formfactor F defined by Eqs. (2, 3) reads 



F(z q ) = - 




In 



z q — ^z q — 1 / 



-i 2 



ITi 



m 



H 



Am 2 



1 



> 1. 



Now we demonstrate how to extract the double logarithmic contribution from the ampli- 
tude using the Sudakov technique. To understand a basic properties of the amplitude in 
DLA we consider first a scalar part of amplitude without numerator and then discuss the 
structure of the numerator separately. The scalar part of the amplitude reads 



F(x q ) 



—2im 2 H 



7T" 



d 4 q 



(ki + q) 2 — m 2 + iO (k 



q) 2 — m 2 + iO 



m 2 + iO 



(5) 



This integral is very simple and can be calculated explicitly. However, we are interested in 
terms containing double logarithms. The most suitable and simplest method of calculation 
was suggested by Sudakov in ||26|| . It consists in a direct integration in the momentum 
space splitting the internal momenta into the parallel and perpendicular components to 
the plane of the external momenta k\ and k 2 . To this end, one introduces new variables 
a and (3 using Sudakov-decomposition for the vector q 



q = aki + (3k 2 + d A q = — —dadf3d 2 q± 



(6) 



We obtain for F (see also consideration of the similar integral in DLA in Ref. [15]) 



F(x q ) = - 

7T 



dad(3 J dp 
o 



x 



where 



\J3 — (x q + p — aj3) + iO] [a + (x q + p — a/3) — zO] [(x q + p — a/3) — zO] 



(7) 



m 



H 



In Eq. ([?D we have performed an integration over the angles in the transverse subspace, 
using the fact that the integrand does not contain a dependence on these variables. It is 
then evident from (0), that the double logarithmic contribution at x q 1 comes from 
the integration region where 

p,x g < \a\,\P\ < 1. (9) 
To carry out the integration over p we may apply the relation 



P- 



[(x q + p — a/3) — iO] x q + p — a/3 



+ in5(x q + p — a/3). 



(10) 
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In the above expression, P means the principal value. 

One can check that the first term in Eq. fllTf ) does not contribute to the DL asymptotics 
of the formfactor and we do not consider it here. We get for the F in DLA 



1 i 



F{x q ) = J *ZL*jLdp9(p)6 (p - (aP - x q )) = J J ^j0(a(3 - x q ). (11) 

— oo —1 —1 

In the last relation, we replace the strong inequality ([]) by usual inequality, limiting the 
integration region in accordance with the logarithmic character of integration. 

Let us now consider the numerator of the amplitude. It is clear from presented deriva- 
tion that only the soft momenta q give contribution to DL. That means that we may 
neglect the dependence on q in the numerator of the integrand: 

Ax r i 

— — d^Sp Y{h + q + m q )(k 2 - q - m q )Y(q + m q ) = 4x q . (12) 

ZftTL j^TTL q 

We see that the numerator is proportional to the second power of the small quark mass. 
One power of the mass appears because of the Yukawa coupling constant and the second 
power is due to the trace. It is important here that only the scalar term of the numerator 
of the t-channel quark propagator (q + m q )/(q 2 — m 2 ) contributes to the amplitude. That 
has a simple physical explanation. The Higgs particle has spin zero and positive parity. 
Therefore, the initial and final total momenta of the subprocess qq —>■ 77 are zero. And, 
the amplitude is suppressed by the mass of the light quark {0{m q / E q )) at high energy 
due to the nonconservation of quark helicity. We see that the quantum numbers of the 
external particles force to have a situation when the t— channel quark propagator behaves 
like a scalar boson propagator, which is singular enough to produce DL. Really , the 
scalar part of the quark propagator is more singular then the vector part in the limit of 
the soft momenta and this fact guarantees the appearance of the new non-Sudakov DL. 
The possibility of the emergency of such DL in the mass suppressed amplitudes due to soft 



quark have been pointed and explained in |J|, where helicity nonconserving amplitude 
77 ~~ ► QQ was studied in DLA. 

Despite the same origin the observed double logarithms in the H — > 77 process demon- 
strate some new peculiarities. For example, this phenomena 1) occurs not in the differ- 
ential distribution or in the special helicity amplitude but rather in the decay width, 2) 
it occurs in the quark loop, but not as a correction to the Born amplitude. Therefore, 
the discussed phenomena is new in a sense that it appears in a very specific and unique 
circumstances, and of course, represents a new manifestation of the non-Sudakov type of 
DL. 

The integration in Eq. (|ll|) is trivial. Taking into account expression fll2"f ) for the 
numerator we get the formfactor F(x q ) in the DLA: 

F{x q ) = Ax q \n 2 {x q ). (13) 

Now we consider a QCD correction to the formfactor F defined by Eq. (|2|). This QCD 
correction was calculated in Refs. [^2, [16|, |24|, |23[ . In the limit of small quark masses the 



correction has a simple form 

1 a s Cp 



F(x q ) = F{x q )(°™- l °°d (l - -^-f- ln 2 (x q ) + ...), (14) 
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where Cp is given by relation 

t a t a = C F I, (15) 

and t a is the color SU(N) group generators in the fundamental representation {Cp = 4/3 
for QCD). There are four different topologies of the two-loop diagrams, which are shown 
in Fig. 1 (b-e). Here we calculate them in the DL approximation by using Sudakov 
technique. 

Let us start from diagram lb and consider first a scalar integral. Here and below we 
always use the Feynman gauge for the gluon fields; the gluon propagator in this gauge is 



k 2 + i0 



(16) 



The scalar part of the F^- lb '(x q ) reads 
F^HxA = 



-2im 2 H 



2 oo 



7T 



(fci + q) 2 - ml + iO {h + Q) 2 - ml + iQ 



x 



(k 2 — q) 2 — m 2 + iO {h,2 — Q) 2 — m 2 + iO q 2 — m 2 + iO (Q — q) 2 — m 2 + iO 
In the Sudakov variables 

q = aki + /3k 2 + q± and Q = Ak\ + Bk 2 + Q_\_, 

we have 

dad(3d q 1 



(17) 



7T 



+ c? -a[3)+ iO] [a + (x q + q 2 - a(3) - iO] 



x 



[(x q + ?-ap)- iO] 



dAdBd 2 Q 



Ti 



B - (x q + Q 2 - AB) + iO 



A 



Q 2 -AB)-iO ((Q-q)2-(A-a)(B-{3) 



iO 



(19) 



where we have switched to dimensionless Euclidian vectors Q and q, noticing that the 
transverse momenta are spacelike. Note that the external integral in Eq. (flSf ) coincides 
with the one loop one. It is then evident that in order to get the DL contribution one 
should integrate over the following region 

f,x q < \a\, |/3| < 1, t<^Q 2 , 

Q 2 ,\a\ < \A\ < 1, Q 2 , < |B| < 1. (20) 

Making corresponding simplifications in Eq. (|i~9f), performing integration over Q and 
g with the help of Eq. (|H]) and taking into account only the absorptive parts of the 
corresponding propagators, we get 



■1 a 



da f d(3 



a J 3 
i -l 



9(a>3 — x 



dA 



dB 



, - (21) 

' >|<|A|<i A J\p\<\B\<i B 6 



6(AB) 



1 



6 



(a) 





Figure 1: Feynman diagrams corresponded to the amplitude (|2]): (a) one-loop diagram 
corresponded to the zeroth order in a s , (b-e) two-loop diagrams corresponded to the first 
order in a s . External wavy lines represent photons, internal wavy lines - gluons, internal 
solid lines - quarks, external solid lines are the Higgs bosons. 
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Because the essential momenta q and Q are soft in DL approximation we can again neglect 
the dependence on them in the numerator in DL approximation. The numerator of the 
integrand then takes the form 



a s Cp Ax q 



Air Am q m 6 H 



Y (ki + q + m q ) 7 A (k 1 + Q + m q ) (k 2 ~Q- m q ) 7 A 



x 



(kz-q- m^j Y (q + m, 
Combining Eqs. (|2l| , |22|) we obtain 



-4x 9 ^. (22) 



F^ b \x q ) = -l^ln^Fix,)^- 1 ^. (23) 

O ATT 

Remaining two- loop diagrams of Fig. l(b-e) do not contribute to the DL asymptotics 
of the amplitude because both the electromagnetic vertex convoluted with a polarization 
vector of real external photon and the quark propagator do not contain double logarithms. 
These properties can be easily checked at the one loop level. Only one large logarithm 
may appear in these structures in the Feynman gauge but not the double one. 

Thus, we have restored the exact result (0). It is clear from (^) that when the mass 
of the Higgs boson is large enough so that in spite of smallness of the QCD coupling 
constant a. 



a s C F . 2 



ln 2 {x q ) ~ 1, (24) 



the perturbative QCD expansion in a s is not valid and we should modify two- loop result. 

In the next section we obtain the quark formfactors F in all orders of QCD in the 
following region of Higgs boson masses 



_«_ lll(lA8) «_ 



«^f^ln(l/i,)«^ifiln 2 (l/i,)~l. (25) 



In this region we expect to get the correct result in the framework of the DLA. 



3 Resummation and Discussion 

The resummed form factor can be presented in the form (see Fig. 2) 

OO A 

F(x ) = 87T1 f 

^ 1' / \ 1/2 / /q_\4 

x tr [S(h + q)V(H -> qq)S{q - k 2 )R^(qq -> 77)] , (26) 

where the amplitude R^ u (qq —>■ 77) is defined in such a way that it contains no two-quark 
cuts in the 77- channel, V(H — > qq) is the exact QCD vertex which contains all possible 
gluon exchanges, S is the exact quark propagator. Double logarithms appear when we 
have consistent logarithmic integrations. Therefore, only region (|S|) can yield the DL 
contribution to the RHS of Eq. (p6|) after integration over q. As it is clear from the above 
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Figure 2: The amplitude in all orders of as- 



discussion we should not consider QCD corrections to the amplitude R fMU (qq — > 77) and 
to the quark propagators S in the DLA, taking them in the Born approximation. One 
has to know only the expression for the H — > qq vertex V(p2,pi) in DLA in the following 
kinematical region 



™> q < \V1l\P2\ < (P2-P1) - 
P2 = h + q, Pl = q-k 2 



m 



(27) 

Let us remind that in general F is proportional to m q . One power of m q appears from 
the Hqq- vertex, the another one - from the trace over Dirac indexes. The physical 
reasons for that were explained in the previous section. There was also noted that one 
of the two powers of m q appears from the quark propagator (q + m q )/(q 2 — m q ) from 
the R^ u {qq — > 77) amplitude. Due to this circumstance we are not forced to consider 
the 0(m q ) corrections to the Hqq- vertex. Therefore, the exact result for this vertex in 
DLA is given by multiplication of the Born expression with the well-known Sudakov form 
factor 



exp 



2vr 



In 



(P2 -Pi) 



I>2 



hi 



Pi 



:28i 



Therefore, the quark form factor F in the DLA can be found from Eqs. (|TTD and (|i~2"l) by 
multiplying the integrand in the RHS of this equation with the Sudakov form factor (|28| ) 
and the result reads 

1 1 



F(x q 



4x n 



-1 -1 



da dp 

w d \ a P 

a p 



Xgjexp 



(o; s C j p/27r) In \a\ In 
The integration over (3 can be easily performed and we have 

— (1 - exp(-Zy(l -y))) , 

y 



(29) 



F x n 



8x q In" 



X n 



(30) 



where Z = a s Cp In 2 (x q )/(2n). We have found that following representation is convenient 



F(x n 



F(x q ) 



(one-loo p) ^^(n+l) ( a s C F 



S r( 2 ^ + 3) 



2tt 



In 2 (x„ 



9 



Fix )( one - lo °p) 



JqJ 6 V 4tt v "V ' 45 V 4tt 



(31) 



The Eqs. (p9|j31| ) are our final results. They correctly describe the quark form factor F 
in the region of the Higgs boson masses where the relations ( p5|) are fulfilled. 

In Fig. 3 we plot relative contributions of two- loop (dashed curve), three- loop (dashed- 
dotted curve), four-loop (dotted curve) order results (0) end exact result in DLA (|30|) 
(solid curve) as a function of Z. All results are normalized on the one-loop result. We 
see that at small enough x q , large Z, the perturbation theory does not work and the 
resummation of DL is mandatory. The middle curve in Fig. 3 corresponds to the exact in 
DLA result. We see that the later slightly decreases the decay width with increasing Z. 

In conclusion, we have discussed a QCD perturbative corrections to the H — > 77 pro- 
cess in the limit of small quark masses in the loops. We have found a new manifestation 
of the non-Sudakov type of double logarithms and have performed one- and two-loop cal- 
culation to find the origin of such double logarithms. Then we have resummed them and 
have derived the exact result for the quark-loop form factor F with double logarithmic 
accuracy ( |29| , |3~T| ) . These results will be important in the situation when the Higss boson 
is heavy enough to guarantee smallness of m q /mH and when the mass m q itself is not 
very small to give visible contribution to the amplitude. Such a situation can be realized 
1) in the Standard Model with Mh ~ 500GeV for the b quark loop and 2) in the Super- 
symmetric extension of the Standard Model, where the coupling of b quark with Higgs 
can be large. For the first application we need to measure the decay width of the H — > 77 
process with very high precision. The feasibility of that depends on experimental side, 
of course, and should be studied separately. We should understand, of course, that at 
very high Higgs masses two-loop electroweak corrections became very important due to 
terms proportional to Gpm 2 H . Fortunately, they have been calculated P5fl . The second 
application seems to be very promising and important phenomenologically. The crucial 
fact here will be the value of Higgs mass and the value of tg(/3). 

The main purpose of this paper was to attract an attention to new non-Sudakov type 
of double logarithms which appears in the H — > 77 process. 

In our opinion, the results presented here may have applications to other cases. One 
of the most interesting is 77 — > qq in the Jz = channel, where such non-Sudakov, 
mass-suppressed DL terms are phenomenologically important, as it was demonstrated in 
Ref. [15]. 
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Figure 3: The relative contribution of the two-loop (dashed curve), the three-loop (dashed- 
dotted curve), the four- loop (dotted curve) and the exact in DLA result for the formfactor 
(j3~T|) (solid curve) as a function of Z . All results are normalized on the one-loop result. 
The one-loop contribution corresponds to one in such a normalization. 
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